Nonequilibrium phase transition in a driven Potts model with friction 
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We consider magnetic friction between two systems of q-state Potts spins which are moving 
along their boundaries with a relative constant velocity v. Due to the interaction between the 
surface spins there is a permanent energy flow and the system is in a steady state which is far from 
equilibrium. The problem is treated analytically in the limit v = 00 (in one dimension, as well 
as in two dimensions for large-g values) and for v and q finite by Monte Carlo simulations in two 
dimensions. Exotic nonequilibrium phase transitions take place, the properties of which depend on 
the type of phase transition in equilibrium. When this latter transition is of first order, a sequence 
of second- and first-order nonequilibrium transitions can be observed when the interaction is varied. 
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I. INTRODUCTION 



Friction is a basic problem in physics with important 
technical implications [1| . Friction between moving bod- 
ies and thus energy dissipation may have magnetic con- 
tributions. This type of phenomenon takes place, for 
example, in magnetic force microscopy when, performing 
a measurement, a magnetic tip is moved over the surface 
of a magnetic material 043 • Recently, magnetic friction 
has been modeled in a simple setup [5|, where two Ising 
models are put in close contact at their surfaces and move 
with a constant relative velocity v. The spins in the sur- 
face layers interact through a coupling Kf and the relax- 
ation process due to phononic and electronic degrees of 
freedom in the material are taken into account via a heat 
bath, at a fixed temperature T, to which all the spins are 
coupled. In this system magnetic friction takes place and 
there is a permanent energy dissipation from the surface 
to the heat bath. The macroscopic motion of the bodies 
represents a permanent perturbation driving the system 
to a steady state which is far from equilibrium. Proper- 
ties of this nonequilibrium state have been investigated 
by Monte Carlo simulations and — in the limit v = 00 
— by analytical methods in different geometries 043 • I n 
the one-dimensional (ID) case there is an order-disorder 
transition at v — 00, but the system stays disordered 
for any finite v in the thermodynamic limit. Fluctuation 
effects introduced by a finite size at v — 00 have been 
recently studied by Hilhorst Q- In the two-dimensional 
(2D) case, when the moving systems are in contact at 
their ID surfaces, the order-disorder transition persists 
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for any v > 0. In all cases the nonequilibrium phase 
transition is found to be of the mean-field type. 

The spin degrees of freedom, which are involved in the 
friction, can have different symmetries and/or different 
types of interactions. As a consequence the equilibrium 
phase transition of the system at v — can be differ- 
ent from that of the Ising model, thereby influencing the 
behavior of the system out of equilibrium, when v > 0. 
In particular, the properties of the nonequilibrium phase 
transition could differ from those of a mean-field tran- 
sition. Therefore, it is of interest to study other driven 
systems with friction and explore the singularities of the 
corresponding nonequilibrium steady states. 

In the present paper we consider a generalization of the 
Ising model with a discrete symmetry, the (/-state Potts 
model [8] , which is equivalent to the Ising model for q = 2 
but displays quite different critical behaviors when q is 
varied. In the equilibrium case there are detailed analyt- 
ical and numerical informations about the physical prop- 
erties of this system. The model has been solved in ID 
for any q, in the large-g limit in 2D, as well as on a fully 
connected lattice. In 2D, exact results are available at 
the phase- transition point 0, which is of second (first) 
order for q < q c (2) — 4 (q > 4). In three dimensions 
(3D), the limiting value g c (3) < 3. Here we investigate 
the g-state Potts model with friction and study the be- 
havior of the steady state as well as the properties of 
the nonequilibrium phase transition for different values 
of q as a function of the velocity v and the friction inter- 
action Kf. In ID, a nonequilibrium phase transition is 
present only at v = 00 where the problem can be solved 
exactly. In 2D we study the quasi-static limit, v — > 0, 
in which case exact results are known about the friction 
force and its singularity at the equilibrium phase tran- 
sition point. At the nonequilibrium case a numerically 
exact treatment is given for v = 00 in the large-g limit. 
These results are compared with numerical simulations 
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replaced by s^. The interaction term is given by 
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FIG. 1: (Color online) Driven interacting Potts models slid- 
ing on each other in (a) ID and (c) 2D, with a translation by 
two lattice constants between interacting sites. For infinite 
relative velocities, the system can be replaced by two decou- 
pled Potts models in (b) ID and (d) 2D, interacting with a set 
of fluctuating variables, attached to one layer of ghost sites 
[green (gray) circles]. These fluctuating variables are equiva- 
lent to an effective surface field with strength hf. 



which are performed at finite v and for different values 
of q. 

The structure of the paper is the following. The model 
is introduced in Sec. UU and the dissipation and the fric- 
tion force in the low v limit are studied in Sec. Mil We 
treat the model in the limit v = oo in Sec. lIVl Numerical 
simulations performed in the 2D case for v and q finite 
are presented in Sec. |V] and the results are discussed in 
Sec. IVII The solution of the self-consistency equation 
for the magnetization in ID is given in Appendix A and 
details about the calculation of the equilibrium surface 
magnetization of the 2D Potts model in the large-g limit 
are given in Appendix B. 



II. MODEL 

We consider two identical Potts models defined by the 
Hamiltonians and H( s \ respectively, which are cou- 
pled at their free surface by a time-dependent interaction 
term V(t). The reduced Hamiltonian of the composite 
system takes the form 



pH{t) = m {s) + m (s ' ] + /3V(t) 



with j3 = l/k B T. 

In ID [Fig. 1(a)] we have 



-K^S(si - s i+ i) . 



(2-1) 



(2.2) 



in terms of the Potts spin variables Si = 0, 1, . . . , q — 1 
and S(n) is the Kronecker delta function. In ftHS s \ Si is 



N 



(2.3) 



where A(t) = vt, and periodic boundary conditions are 
used, i + N = i. 

In 2D [Fig. 1(c)] the Potts variable s(i,j) are attached 
to the sites of a square lattice. The Hamiltonian of 
the Potts model is then 

JV L 

m (s) = -KiY.Y, 6 ^*-*** 

«=1 3=1 
N L-l 

«=1 3=1 

and similarly for ^HS S >. The couplings in the horizontal 
and in the vertical directions, K\ and K2, can be different 
and the interaction term takes the form: 



N 



\+A(t) 



(2.5) 



We set periodic boundary conditions in the horizontal di- 
rection, 8i+N,j = Si j as well as in the vertical direction, 
which means that there are two equivalent sliding inter- 
faces in the system. The problem is generally studied 
in the thermodynamic limit. Note that the 2D problem 
with L = 1 (i.e., with one layer in each subsystem) is for- 
mally equivalent to the ID problem with s(i) = s(i, 1). 



III. DISSIPATION AND THE FRICTION 
FORCE FOR v -> 

To see the relation of the above model with mag- 
netic friction we follow Ref. [5J], couple the system to 
a heat bath of constant temperature, T, and study its 
nonequilibrium properties by Monte Carlo simulations. 
In our case the relaxation kinetics is governed by the 
heat-bath algorithm Measuring the energy differ- 

ence, AE = E' — E, between the original (E) and the 
flipped (E 1 ) configurations it is found that AE < ( 
i.e. energy is dissipated to the heat bath, which shows 
the presence of magnetic friction in our system). More 
detailed investigations are performed for the 2D prob- 
lem with L x L spins (i.e., for N = L with L x L/2 
subsystems) and with periodic boundary conditions, in 
which case the dissipated energy per spin during time 
t, Ai?bath(i) / L 2 , has been measured for different rela- 
tive velocities, v, and temperatures, T. We illustrate the 
time dependence of AEb at h(t)/L 2 in Fig. |for d = 1 
and for Kf = K\ = K2 = 1/T at a temperature 
T = l.lT c (g), above the bulk phase transition point, 
T c (q) = l/[ln(l + y/q)], for different values of q. As for the 
Ising model 5] with q = 2 the dissipated energy grows lin- 
early in time, Ai?bath(0 = Pt, and the dissipation rate, 



3 




400 



0.25 



0.20 



0.15 



0.10 



0.05 



0.00 



o. 




q=2 

□ — □ ^=5 



0.6 0.8 1.0 1.2 1.4 1.6 1.8 
T/TJq) 



FIG. 2: Accumulated energy per spin, which is dissipated by 
the heat bath during time t in a 80 x 80 lattice for different 
values of q. The sliding velocity is v = 1 and the temperature 
is T = l.lT c (g), see the text. 



FIG. 3: The magnetic frictional shear stress, F/L, as a func- 
tion of temperature in a 80 x 80 lattice for different values of 
q and a sliding velocity v = 1. In the thermodynamic limit 
for q = 2 and 3 the derivative of F/L is singular at T c (q), 
whereas for q — 9 and 16, F/L has a discontinuity, which is 
also indicated in the figure. 



P, depends linearly on the velocity for small v. P = Fv, 
where F is the friction force. F is found proportional 
to the length of the cut, L, so that the magnetic fric- 
tional shear stress F/L can be used to characterize the 
magnetic friction. 

We have measured the magnetic frictional shear stress 
at different temperatures. Data are presented in Fig. [3] 
for q = 2, 3, 9 and 16. In a finite system with 80 x 80 spins 
F / L shows some kind of extremal behavior in the vicin- 
ity of the bulk phase-transition temperature, T — T c (q). 
For q = 2 and 3, having a second-order equilibrium tran- 
sition, the derivative of F/L with respect to T shows 
a maximum, whereas for q — 9 and q = 16, having a 
first-order equilibrium transition, F/L itself is maximal 
around that point. In reality one is interested in the be- 
havior of the system in the thermodynamical limit, in 
which case we define / = lim£_ i . 00 F/L. In the limit 
v — > 0, the sliding velocity is so slow that the system has 
time to relax to equilibrium between successive relative 
moves of the two subsystems. Then one can express / 
as the difference between two equilibrium spin-spin cor- 
relation functions [5|. In the original configuration E/L 
is proportional to the nearest-neighbor correlation func- 
tion, C\ — (<5(si,i — s'ii)), and after the displacement 
E' / L is proportional to the next-nearest-neighbor corre- 
lation function, C% = (S(si t i — s\ 2 )}, so that the magnetic 
frictional shear stress is given by f — (G% — Ci)/2, where 
the division by 2 is due to the two equivalent sliding sur- 
faces. For the Ising model these correlations are known 
and thus f(T) can be calculated exactly. Its deriva- 



tive at T = T c shows a logarithmic singularity: 



dT 



^d-V2)ln 



l-L 



(3.1) 



Comparing this result with the numerical findings in 
Fig. [3] we can say that in the thermodynamic limit the 
slope of the curve at T = T c and q = 2 is divergent for 
v —> 0. We expect that some kind of divergency will stay 
for a finite v also and the finite slope in Fig. [3] is a finite- 
size effect. One can generalize this result for q = 3 and 
more generally for systems having a second-order equi- 
librium transition. Since the near-neighbor correlations 
have the same type of singularity as the energy density, 
we obtain in the v — > limit: 



dT 



1 - 



T 



(3.2) 



where a is the specific heat critical exponent of the sys- 
tem (for the q = 3 Potts model it is a = 1/3). On the 
contrary, for q — 9 and q — 16 and more generally for sys- 
tems with a first-order equilibrium transition, the energy 
density, and thus the magnetic frictional shear stress has 
a discontinuity at T = T c . Consequently the true behav- 
ior in Fig. |3j for infinite systems is a jump for q = 9 and 
16, at least in the slow displacement (v — > 0) limit. This 
limiting behavior is indicated by the dashed and vertical 
lines in Fig. [3] 

In the experimentally relevant situation the sliding ve- 
locity is finite and the system under investigation is out 
of equilibrium. In this already demonstrated for 
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the Ising model [6j, a nonequilibrium phase transition 
takes place in the system. The nonanalytical behavior of 
the friction force in Fig. 3 strongly suggests the existence 
of a nonequilibrium phase transition for q > 2 as well. In 
the following we study this nonequilibrium phase transi- 
tion, first in the infinite velocity limit, in which case the 
mean-field treatment is exact, and afterwards for finite v 
by Monte Carlo (MC) simulations. 



IV. SOLUTION AT INFINITE RELATIVE 
VELOCITIES 

A. Potts model with a fluctuating variable 

The solution of the problem becomes simple at v = oo 
as noticed already for the Ising model in Ref. @ . Here we 
use a generalization of the same argument for the q-state 
Potts model. The basic observation is that at v = oo, the 
spin variables s^i and s' i+A ^ 1 in the interaction term 
of Eqs. (I2.3[) and (|2.5p are uncorrelated. Consequently, 
for the subsystem (s), the spin s' i+A , t ) can be replaced 
by a randomly chosen spin from the surface layer of the 
subsystem (s'). The same effect on Sj ( x is obtained if 
the surface spin is coupled to a fluctuating Potts variable 
/ij = 0, 1, . . . , q — 1, with the following constraints on the 
mean- values 



(5( fH -a))={6(s'-a)) = 



[S(a)q — l]n%f + I 



(4.1) 



where m/ is the mean value of the magnetization in the 
two interface layers and a = 0,1, . . . ,q — 1 is a Potts 
variable. Thus the interaction term in Eq. (|2.5|) can be 
replaced by the following effective interaction for the sub- 
system (s) (see Fig. [T]) 

S(s h i - s-+A(t),i) S ( s iA - Mi) , ( 4 -2) 
with the probability distribution 
1-m/ 



S(/j,i)mf 



(4.3) 



for the fluctuating variable. 

In the next step we integrate out the fluctuating vari- 
able Hi and replace its effect by an external field acting 
on Si i. 



B. Effective external field 

Here we consider a Potts model on a general lattice, 
with Hamiltonian T-Lq involving the Potts variables s k - 
These variables are interacting with a set of fluctuating 
Potts variables fx k , through the couplings K k . The total 
reduced Hamiltonian is given by 



fiW-n = PH-a - ^2 K k S(s k - n k ) 



(4.4) 



where the /i k are distributed as in Eq. (|4.3[) with a mean 
value, n%k- 

The fluctuating variables fi k can be traced out from 
the partition function Z^. Making use of the identity 

e K k S( Sk -„ k ) = 1 + {e K k _ 1)s{sk _ Mfe) ^ (4 5) 

one obtains 

= Tr s ^e- 0n " = Tr 5 e-^°Tr Al JJe Kfe ' 5 ^-^ 

= ^Qc(X /c ,m fe ) 

k 

xTr s e- fSn °l[[l + d(K k ,m k )5(s k )] , (4.6) 



where: 



c(K kl m k ) = 
d(K k ,m k ) = 



1 + {e Kk - 1 



1 - m k 



q(e Kk - l)m k 



q+{e K * -l)(l- mfc ) 



(4.7) 



The expression for the partition function in Eq. (|4.6[) 
can be compared to that of a Potts model, coupled to a 
static magnetic field h k , and defined by the Hamiltonian: 



(3H cq = (3n Q -J2 hkS ^ ■ 



(4.8) 



For this model the partition function reads: 

Z cq = Tr s e- m ^ =Tr s e~ m °Y[e hkS ^ 

k 

= Tr s e-? n ° {J [1 + (e hk - l)6(s k )] . (4.9) 
fc 

If we fix the values of the field h k to h k such that e hk — 1 = 
d(K k ,m k ), that is 



h k = In il 



q(e Kk - l)m k 



q+{e K * -l)(l-m fc ) 



(4.10) 



then the partition function in Eq. (I4.C5[) can be expressed 
with the equilibrium partition function as: 



Zn = Y\_c(K k ,m k ) x Z eq (h) 



(4.11) 



C. Application to the driven Potts model 

Now we turn back to our original problem where two 
interacting Potts models are moving with a constant rela- 
tive velocity v and have a magnetization mj at the inter- 
face in the stationary state. In the limit v = oo this sys- 
tem can be replaced by two noninteracting Potts models 
interacting with a set of fluctuating variables, attached 
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to one layer of ghost sites as shown in Fig. [T] Integrat- 
ing out the degrees of freedom associated with the ghost 
sites, the effect of one subsystem on the other is equiv- 
alent to a static effective surface field with strength hf. 
Using the formalism of Sec . IIVB1 we have = m/ and 
Kk = Kf for all surface sites k, whereas Ho is the Hamil- 
tonian of the subsystem defined in Eqs. (|2.2p and (|2.4I) . 
The effective surface field is also the same for all surface 
sites, hk — hf, and follows from Eq. (|4.10l) so that 

e h f = l + n(m f ,K f ) 

= f^L, (4-12) 



-{K f )m f 



with 



r(K f ) 



e K f-l 
e K f + q - 1 



(4.13) 



The surface magnetization in the equilibrium system is 
obtained as 



m/, eq 



q d\nZ cq 
N dh f 
9-1 



and must satisfy the self-consistency equation 
TO /,eq [hf(K f ,m f )] =m f , 



(4.14) 



(4.15) 



with the appropriate values of the subsystem interac- 
tions. 

When the transition is of second order, the transition 
point satisfies the condition dmf iCq /dmf\ mf= o = 10- 
Using Eqs. (|4.12|) and f)4. 15[) we arrive at 



X { f °l q \ c qr(K fc ) = l, 



(4.16) 



where X/°l q = d m f,e<i/dhf\h f =o is the zero-field surface 
susceptibility of the g-state Potts model. 



D. Analytical solution in ID 

The magnetization of the ID Potts model in the pres- 
ence of an external field hf can be calculated by the 
transfer-matrix method. This is explained in Appendix lAl 
where, using the relation of Eq. (|4.12[) between the mag- 
netization rrif and the effective field hf, a closed form 
for the self-consistency condition in Eq. (|4.15|) is derived. 
Besides the trivial solution 



mf =0, 



(4.17) 



the other solutions are given by two roots of the cubic 
equation 



a 3 m f 



a 2 m 2 f 



airnf + cto = , 



(4.18) 



with coefficients: 

a = _(g_2)[2r(e JC -l)+g(r-l)] 

ai = -2r(q - 2) 2 (e K - l) - q(q - 1) 

+ r 2 { [q(e K -l)+ 4( 9 - 1)] (e K - l) + q(q - 1)} 

a 2 = T { q -2){e K -l) 

x {T[q(e K -l)+ 4( g -l)] +2( 9 -l)} 

a 3 = -(q - l)r 2 (e K - l) [q (e K - l) + 4(q - 1)](.4.19) 

As explained in Appendix El the third root does not sat- 
isfy the self-consistency condition in Eq. (|A5j) . We have 
discarded this nonphysical root after a direct substitution 
into Eq. (|A5]) . 

The structure of the cubic polynomial is different for 
q = 2 (Ising model) and for q > 2. For the Ising model, 
due to the up-down symmetry, the even coefficients in 
Eq. (|4.18p are vanishing, ag = a 2 = 0, and the nontrivial 
solutions are: 



(±) , / r 2 e 2 *-l 



(4.20) 



We have checked that these solutions are the stable ones 
below the critical temperature T c which is given by the 
condition: 



-(K fc ) 



= 1 



(4.21) 



Note that for the Ising model r(Kf) = tanh(if//2V, thus 
we recover the result previously obtained in Ref. [6j] . As 
the critical point is approached the nonequilibrium mag- 
netization is vanishing continuously with a critical ex- 
ponent (5(q = 2) = 1/2. The temperature dependence 
of the nonequilibrium magnetization for the symmetric 
case, Kf = K, is shown in Fig. |U Here we use the tem- 
perature parameter & q , defined as: 



V2{e K 



1) 



(4.22) 



For the symmetric Ising model at the critical point we 
have 6 2c = 1 [see Eq. (|43T]) ]. 

Now we turn to the solution of the non-Ising case, 
q > 2, looking for those roots of the cubic polynomial 
in Eq. (|4.18[) which are the nontrivial solutions of the 
self-consistency equation. We know from Cardano's for- 
mula that the structure of the real solutions of a cu- 
bic polynomial depends on a discriminant defined as 
Disc = Q 2 + P 3 , with: 



P 



1 

27 



«2 
«3 

«2 
«3 



ai 
3a 3 

aia 2 



6a 2 



2a 3 



(4.23) 



For Disc > 0, which happens when Q q > Q qc , the poly- 
nomial has one real root and there is no nontrivial so- 
lution of the self-consistency equation. On the contrary, 
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K=K/2 
K f =K 

K r 2K 

q=6 
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FIG. 4: (Color online) Temperature dependence of the mag- 
netization m/ of the g-state Potts chain with friction in the 
limit v = oo for the symmetric model with K = Kf. The sta- 
ble (unstable) solutions of the self-consistency equation are 
denoted by full (broken) lines. The transition is of second 
order for q = 2 (Ising model) and first order for q > 2. In 
the latter case there is a hysteresis: on heating (cooling) the 
transition is at Q qc (Qqi = 1). The inset gives the deviation 
of the transition point O qc from the Ising value &2c = 1, as 
well as the value of the jump in the magnetization on heating, 
m/ c , as a function of q. 

Disc < in the low temperature region such that < 
Q q < Q qc and the polynomial has three real roots. For 

the two nontrivial solutions, such that > we 

have checked that mi > in the whole region whereas 

m { - ] < (m { f - ] > 0) for 9, < @ ql (Q ql < Q q < Q qc ). 

Thus Qqi is defined by the condition m\ ' = 0. The two 
nontrivial solutions annihilate at the transition point Q qc - 

We have studied the stability of the solutions as well 
as their domains of attraction by considering the self- 
consistency Eq. (|A5[) and varying rrif on its right-hand 

side. The solution ^ is always unstable, on the con- 
trary is always stable. The region of attraction of 
the latter solution isl>TO/>0(l>m/> m\ ') for 

Q q < Qqi (Qqi < Q q < Q qc )- The trivial root m ( p = 
is unstable for Q q < 9 ?1 and stable for Q q > Q q i with 

the region of attraction m/ < (0 < m/ < 1) for 

6,1 < Qq < Qqc (Qq > Qqc)- 

The magnetization of the symmetric model for differ- 
ent values of q is shown in Fig. [4] as a function of Q q . For 
q > 2 the nonequilibrium phase transition is first order 
and the magnetization shows a hysteresis. Starting from 
the ordered phase with Q q < Q q \ and heating the system, 



FIG. 5: (Color online) Influence of the coupling Kf between 
the two chains on the temperature dependence of m/ for q = 
6. The stability of the ordered phase increases with Kf (from 
left to right). 

the magnetization remains given by the nontrivial stable 
solution rrS^ > until Q qc where it jumps to to^°' = 0, 
the trivial solution. In the reverse process, cooling the 
system from the disordered phase with Q q > Q qc , the 
magnetization remains vanishing (m*p = 0) until Q q \ 

where it jumps to the nontrivial solution > 0. The 
location of the transition point Q qc and the value of the 
magnetization jump mf c — rnf(Q qc \_) ~ m/(6 gc | + ) are 
shown in the inset of Fig. |4] for the symmetric model as 
a function of q > 2. 

For q > 2, gl corresponds to a vanishing nontrivial 

solution \ thus to ao = according to Eq. (|4.18|) . 
Using Eqs. (|4T3|) and (|422|) to express a in Eq. (|4T9)) 
leads to: 

It follows that in the symmetric case, Kf/K = 1, one 
obtains 8 gl = 1 as shown in Fig. @] When Kf/K in- 
creases, the ordered phase becomes more stable and Q q i 
increases, too. 

The influence of the coupling ratio Kf/K on the tem- 
perature dependence of the nonequilibrium magnetiza- 
tion is shown in Fig. [5] for q = 6. 

We have the following asymptotics for the symmetric 
model. When q is close to 2, the critical point and the 
magnetization discontinuity behave as: 

Q qc ^ l + i(g-2) 2 , m /c ~<7-2. (4.25) 
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When O g < O gc , the magnetization approaches the lim- 
iting value nifc with a square-root singularity: 



m f - 171 fc ~ V ®<?c - 6 9 • (4.26) 
For large-g values the transition point is located at 

1 - aq- 1 ' 6 







l + ag- 1/6 , 

(4.27) 



and the magnetization discontinuity is given by 

m /c w l-&g- 1/6 , (4.28) 
with a « 0.9447 and 6 « 0.629, thus a/6 w 3/2. 

E. Solution in 2D for large-g values 




The key point of the solution of the driven system in 
the limit v = oo is the knowledge of the equilibrium sur- 
face magnetization as a function of temperature and sur- 
face field. For the Potts model in 2D, analytic results 
about the surface magnetization are known for q = 2 (Tlj 
(Ising model) as well as in the large-g limit [r2Hl4| . The 
analysis for the Ising model has been performed in Ref. [|| 
and here we consider the Potts model in the large-g limit. 

The 2D Potts model is defined in Eq. ([2~4]l with the 
time-dependent interaction term given in Eq. (|2.5p . We 



treat the problem in the strongly anisotropic limit 15 1 
where the horizontal coupling K\ — > oo, the vertical cou- 
pling if 2 — > 0, while the ratio J = K-z/K* remains con- 
stant. Here K\ is the dual coupling defined through: 



1 



,A7 



1 



(4.29) 



Then the column-to-column transfer matrix of the non- 
interacting system in Eq. (|2.4p takes the form T = 
cxp(— K*H), where H is the quantum Hamiltonian [ll|: 



L-l 



L q-l 



(4.30) 



In the last term, M, 



j=i fe=i 



3.j is a spin-flip operator such that 
oj — a j t fc,mod(g). The quantum Potts model de- 
fined in Eq. (|4.30|) has a quantum phase transition at 
J = J c = 1 in the thermodynamics limit. This tran- 
sition is of second order for q < 4 and first order for 
q > 4. The interaction between the two driven systems 
in the limit v — oo is represented by a static surface 
field, the strength of which is obtained from Eq. (|4.12j) 
as hf — Kfnif in the strongly anisotropic limit, Kt — > 0. 
Thus, using the parametrization Kf = K^K/^/q, the 
Hamilton operator in Eq. (I4.30P is supplemented by a 
surface-field term: 



V 



-hS(si 



K* 



(4.31) 



FIG. 6: (Color online) Temperature dependence of the in- 
terface magnetization of the 2D Potts model with friction in 
the v = oo and large-g limits, using the strongly anisotropic 
(Hamiltonian) version of the problem, for values of the in- 
terface coupling k increasing from left to right. When the 
strength of the interface coupling increases, there is a critical 
value k c above which the system becomes bistable with an 
unstable solution (broken line) between two stable ones (full 
lines). The interface magnetization always vanishes in the 
bulk disordered phase, 6 > 0. 



The surface critical behavior of the quantum Potts model 
in the large-g limit has been studied in Refs. [12H14| and 
the results are summarized in Appendix [B] In the large-g 
limit, corrections of the order of g -1 / 2 are taken into ac- 
count and the distance from the critical point is defined 
as J c — J = 6/-Jq, where 9 plays the role of a reduced tem- 
perature. Using these results we can calculate the equi- 
librium surface magnetization m/ )eq [#, hf(n, mf)] and 
the nonequilibrium magnetization is deduced from the 
self-consistency condition in Eq. f|4. 15|) . We have solved 
the self-consistency equation for the interface magneti- 
zation numerically exactly (l~7| . The dependence of m/ 
on the reduced temperature 8 for different values of the 
interface coupling k is shown in Fig. [5] 

For k = we recover the magnetization at a free sur- 
face, which vanishes linearly at the bulk critical temper- 
ature C — 0. With increasing interaction k the interface 
magnetization increases for 9 < 0, but it is always van- 
ishing in the bulk disordered phase, 9 > 0. This behavior 
is due to the fact that the equilibrium surface magneti- 
zation of the model mi — above 9 C , whatever the value 
of the surface field, as seen in Appendix |B] This is a 
peculiarity of the system in the large-g limit which will 
be discussed further in Sec. IV CI 

In the ordered phase, the shape of the magnetiza- 
tion curve qualitatively changes with increasing k. At 
a critical value k' c — 2.104, its slope diverges when 



8 



9 — 9' c = —0.13037. In the vicinity of this critical temper- 
ature the nonequilibrium interface magnetization shows 
a power-law singularity m(«4, 9) — m{n' c , 9' c ) ~ \9— 9' c \ 1 / 3 . 
Increasing the strength of the interface coupling further, 
k > n' cl a bistability occurs. Here in a finite range of 
temperature, 9% < 9 < 62 < 0, there are three solu- 
tions of the self-consistency equation, which are denoted 

(1) (2) (3) (2) 

by rrij < rrif < rrif . Among these rrif is unstable, 

whereas is stable for perturbations in the region 

to/ < nrSp and rrSp is stable for to/ > to^ 2 -*. Conse- 
quently starting at sufficiently low temperature, 9 < 9\ 
and heating the system, its interface magnetization will 

stay on the upper part of the curve and follows the solu- 

(3) 

tion rrif in the range 9\ < 9 < 92 . At 9 = 6*2 it jumps to 
the lower part of the curve, which is the continuation of 
the solution to^ 1 " 1 . This jump being finite the transition 
is of the first order. When #2 < 0, heating the system 
further the interface magnetization vanishes at 9 C = 
linearly. In the reverse process we start in the disordered 
phase, 9 > 0, and cool down the system. Then the inter- 
face magnetization increases linearly below 9 C = and 
follows the solution mf , in the range 62 > 9 > 9\. At 

9 = 9\ it jumps to the upper part of the curve, which 

(3) 

is the continuation of the solution m\ . Consequently 
there is a hysteresis in the temperature dependence of 
the nonequilibrium magnetization. 

The limiting value 92 first increases with k until k = 
[27 

k* = W — R3 2.598 above which it stays at #2 = 0. For 
k > k* on heating the nonequilibrium interface mag- 
netization jumps directly from rrif (0) > to to/ = 
at 9 = [18|. Thus, on heating, for strong interaction 
k > k* , there is a first-order nonequilibrium transition 
at 9 = 0. In the reverse process, on cooling, the second- 
order transition at 9 — is followed by a first-order one 
at 9 = 6>i. 



V. MONTE CARLO SIMULATIONS OF THE 2D 
SYSTEM 

In the present section the results obtained previously in 
the limit v = 00 are confronted with Monte Carlo simula- 
tions. In particular, we want to see how the fluctuations, 
introduced by a finite relative velocity between the two 
driven systems, influence the properties of the nonequi- 
librium interface magnetization and its singular behavior 
at the phase-transition point. For the Ising model with 
q = 2 these questions have been studied in Ref. [gj]. In 
ID the ordered phase is suppressed for any finite value of 
v. On the contrary in 2D the nonequilibrium fluctuations 
introduced by the finite velocity are found to be irrele- 
vant and the nonequilibrium phase transition is described 
by the same (mean-field) critical exponents. 

Concerning the Potts model in ID the nonequilibrium 
fluctuations should destroy the ordered phase for any fi- 



nite value of g, thus one does not expect any nonequi- 
librium phase transition for v < 00. However, the ques- 
tion is more delicate in 2D, where the properties of the 
nonequilibrium phase transitions in the limit v — 00 are 
different for q = 2 and for large-g values. These changes 
are likely to be related to the surface-field dependence 
of the surface magnetization in the equilibrium systems 
[see the self-consistency condition in Eq. f|4. 1 5[) ] . which 
is known to depend on the order of the bulk equilibrium 
phase transition. Therefore it is advisable to study the 
two regimes in the Monte Carlo simulations. Thus we 
treat successively the case q = 3, where the equilibrium 
transition is second order, in Sec. lVBl and the case q = 9, 
where it is strongly first order, in Sec. IV CI 

A. Method of simulation 

We simulate systems consisting of L x L spins (i.e., for 
N = L with L x L/2 subsystems), with L ranging from 
80 to 320. We use symmetric couplings between pairs 
of spins, with K\ — K2 — l/T, but allow for a varying 
coupling strength Kf across the cut separating our two 
Potts systems, with Kf ranging from 0.25/T to 2/T. In 
the following, we denote the ratio Kf /Kj, = KfT as n. In 
the horizontal direction (i.e., the direction parallel to the 
interface) periodic boundary conditions are used. In the 
vertical direction we use both periodic boundary condi- 
tions (yielding two interfaces) and open boundary condi- 
tions and find that, in general, the local quantities close 
to the interface are independent of the vertical boundary 
condition. Small deviations, if any, only show up for the 
smallest system size. The data discussed in the following 
have been obtained with periodic boundary conditions in 
both directions. 

Our main focus is on the magnetization profile and, es- 
pecially, on the magnetization close to the interface sep- 
arating the two g-state Potts systems that move with the 
relative constant velocity v. The magnetization of row i 
is given by 

(m(ty = (qN m /L-l)/(q-l) (5.1) 

where N m = max(Ao, Ni, N q ). Here N q is the aver- 
age number of spins in state q in row i. Obviously, we 
have that to/ = (m(l)) = (m(L/2)) = (m(L/2 + 1)) = 
(ML))- 

For the Monte Carlo updates we use the standard heat- 
bath algorithm as our single-spin flip algorithm. For the 
corresponding 2D Ising model [6] it was shown that qual- 
itatively the results obtained in simulations are indepen- 
dent of the update scheme. This is different when looking 
at specific quantities, as, for example, the phase transi- 
tion temperature, which do depend on the chosen algo- 
rithm. In the following we restrict ourselves to a quali- 
tative discussion of the properties of our nonequilibrium 
system. 

We implement the sliding of one half of the system 
with respect to the other half in the same way as in 




FIG. 7: (Color online) Interface magnetization for the q — 3 
model, with (a) v = 10 and various values of k and (b) k = 2 
and various values of u. In all cases a continuous surface 
phase transition is observed. The data shown here have been 
obtained for a system containing 160 x 160 spins. 

Refs. [5, 6]. When simulating a system with sliding ve- 
locity v we translate the upper half of our system by 
one lattice constant after L 2 jv random sequential single 
spin updates. One Monte Carlo step therefore consists of 
L 2 single spin flips and v translations. In our numerical 
study, we varied v between 1 and 80. 

Based on the results discussed in the previous sections, 
we expect to observe a discontinuous change of the mag- 
netization close to the cut for q > 4. To observe a pos- 
sible hysteresis we made both heating and cooling runs. 
Starting the heating (cooling) runs with a fully ordered 
(disordered) initial state, we typically let the system relax 
for 2 10 5 Monte Carlo steps at the first temperature be- 
fore starting the measurement. After averaging over typ- 
ically 10 5 steps, we changed the temperature and let the 
system relax for a few 10 000 time steps before starting 
the measurement at the new temperature. We carefully 
monitored our system to detect the possible presence of 
a hysteresis. This procedure was repeated at least ten 
times and the data discussed in this Section result from 
averaging over these independent runs. 



B. Results for q = 3 

The equilibrium q = 3 Potts bulk system exhibits a 
continuous phase transition at the temperature T c (q = 
3) = l/ln(l + v3) « 0.995, similar to the Ising model 
which can be viewed as the q = 2 Potts model. We 
show in Fig. [7] the temperature dependence of the inter- 
face magnetization for various coupling strengths k and 
various sliding velocities v. In all cases we observe a 
continuous boundary phase transition at a critical tem- 



FIG. 8: (Color online) Interface magnetization for the q = 9 
model, with (a) v = 10 and various values of k and (b) k = 1 
and various values of v. The order of the surface transition 
changes from continuous for weak interface couplings to dis- 
continuous for strong interface couplings. Note the changes of 
temperature scales. The discontinuous character of the phase 
transition is revealed by the presence of a thermal hysteresis 
[filled (open) symbols on heating (cooling)]. 

perature that depends on both v and n. Thus an in- 
crease of the coupling strength across the cut yields an 
increasing strength of the effective surface field that sta- 
bilizes the interface magnetization against thermal fluc- 
tuations. An increase of the boundary phase transition 
temperature is also observed when increasing the slid- 
ing velocity, which is similar to what is observed for the 
Ising model @. In fact, the q = 3 Potts model be- 
haves in every aspect like the Ising model. Both mod- 
els have a continuous bulk phase transition, in both 
models the boundary phase transition is continuous and 
of mean-field type [we checked the mean-field charac- 
ter of our transitions by studying the effective exponent 
/3 e g = din m/ /d ln(T c — T) and found that this exponent 
tends to 1/2 when approaching T c ]. 

C. Results for q — 9 

To see whether the scenario obtained for the v = oo 
and large-q limits is generic for values of q > 4, we have 
studied intensively the case q = 9. For that value of q 
the bulk system undergoes a strong first-order transition 
at the temperature T c (q = 9) = l/ln(l + Vd) w 0.721. 

Our main findings are summarized in Fig. [5] Fixing v 
and changing the interface strength k reveals two inter- 
esting features, see Fig.|5[a). First we note that for small 
values of k the boundary phase transition is continuous 
and takes place at a temperature that is comparable to 
the temperature of the bulk transition. For larger val- 
ues of k, however, the boundary transition is discontin- 
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uous as revealed by a thermal hysteresis. We therefore 
have the interesting situation that a continuous and a 
discontinuous surface transitions are separated by a tri- 
critical point. This behavior is in full agreement with 
the v = oo and large-q scenario discussed in Sec. HID. 
There is, however, also a remarkable difference between 
the q = 9 system and the large-*? case. Whereas in the 
large-g case the interface magnetization is strictly zero 
above the equilibrium phase transition temperature, ir- 
respective of the value of the interface coupling, for q = 9 
we hnd for large values of n a hnite surface magnetiza- 
tion above T c (q = 9), where the bulk is disordered, with a 
subsequent discontinuous surface transition at some tem- 
perature T s > T c (q — 9). This unexpected behavior 
can be understood by mentioning that the correspond- 
ing equilibrium semi-inhnite system has a finite surface 
magnetization in strong enough fields for temperatures 
above T c (q = 9) 0. 

Interestingly, the surface transition temperature for a 
fixed value of n is largely independent of the magnitude 
of the sliding velocity v. This is shown in Fig. |8jb) for 
k = 1. Obviously, the limit v = oo is approached very 
rapidly and already modest values of v yield results that 
are very close to those expected for large values. 



VI. DISCUSSION 

In this paper we have studied the magnetic contribu- 
tion to the friction in the g-state Potts model in which 
two interacting systems are moving with a constant rela- 
tive velocity v. During the movement of the macroscopic 
bodies there is a permanent energy flow into the heat 
bath, resulting from the friction force. The system is 
thus driven into a nonequilibrium steady state, which 
can show order-disorder phase transitions as the temper- 
ature, the strength of the interaction or the velocity are 
varied. The Potts model, with its rich critical behavior 
at equilibrium, depending on the value of q, is a suitable 
system to study under nonequilibrium conditions. 

We have studied the phase diagram and the phase tran- 
sitions in this nonequilibrium system for different values 
of q in ID and 2D by analytical and numerical methods. 
In ID long-range order is present only in the v = oo limit, 
where fluctuations are completely suppressed, and the 
problem is solved exactly using the mean-field method. 
For finite v in a system of finite extent, N < oo, one 
expects cross-over phenomena in analogy to the q = 2 
case @ . The nonequilibrium phase transition in this sys- 
tem is of second order for the Ising model (q = 2) but 
of first order for q > 2. In the latter case there is a 
hysteresis: the jump in the nonequilibrium magnetiza- 
tion takes place at different temperatures on heating and 
cooling the system, respectively. This type of behavior 
is expected to take place for other driven ID models for 
v = oo also, provided the up-down symmetry of the local 
order parameter is absent. 

In 2D the nonequilibrium phase diagram is found to be 



more interesting and more exotic. Here the quasi-static 
limit, v — > 0, is different from the true nonequilibrium 
case v > 0. In the quasi-static limit the friction force 
has the same type of singularity at the equilibrium phase 
transition point as the equilibrium energy density. On 
the contrary, for any finite v > the singularity in the 
steady state at the phase transition point is controlled 
by another fixed point, in which the critical exponents 
are mean-field like. More detailed results are obtained 
in the large-g limit, where for large-w the problem is 
solved exactly using the mean-field method. Here the 
phase transition is of second order for sufficiently weak 
interaction. With increasing interaction a second phase 
transition takes place, at a lower temperature, which is 
continuous for a critical value of the coupling and dis- 
continuous for larger couplings. This latter transition is 
accompanied by a hysteresis. Numerical studies of the 
q = 9 state Potts model with finite velocity have shown 
a similar scenario as described above for large q. 

On the contrary, numerical results obtained for the 
q = 3 model have shown just one continuous transition 
for any value of the couplings, which is of the mean-field 
type. This is similar to the scenario found for the Ising 
model. This difference in the phase diagram is expected 
to be related to the nature of the corresponding equilib- 
rium phase transition. If the equilibrium phase transi- 
tion is continuous, which happens for q < 4, the surface 
phase transition is continuous also and the zero-field sur- 
face susceptibility is nonzero even above the bulk transi- 
tion temperature. In this case, a second-order nonequi- 
librium transition is expected to take place at a tempera- 
ture which is higher than the equilibrium bulk transition 
temperature. On the contrary, if the equilibrium phase 
transition is first order, which is the case for q > 4, the 
surface transition is usually second order, a phenomenon 
which is known as surface-induced disorder [19h23| |. In 
this case the zero-field surface susceptibility is zero at 
and above the bulk transition temperature and one needs 
a finite surface field, h s > h sc (q) > 0, to have a nonvan- 
ishing surface magnetization at and above the transition 
point [24] . This noncontinuous surface field dependence 
of the surface magnetization is responsible for the differ- 
ent scenario in the nonequilibrium system for q > 4, in 
particular, for the first-order transition. The hysteresis, 
which accompanies the first-order transition, is due to the 
nonequilibrium nature of the process and the area of the 
hysteresis loop is proportional to the energy dissipated 
during the transition. 

Since the properties of surface- induced disorder are ex- 
pected to have the same type of discontinuous surface 
field dependence for any equilibrium first-order transi- 
tions in 2D and 3D [20l | , the same type of nonequilibrium 
scenario, which we have found for the 2D Potts model 
with q > 4, is expected to take place in these systems. 

We close our paper with some remarks about the possi- 
ble occurrence of nonequilibrium phases and phase tran- 
sitions in other driven systems. As we have seen in our 
study, there is an intimate connection between the equi- 
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librium surface critical behavior of these systems and the 
nonequilibrium states with friction. This relation is ev- 
ident in the limit v = oo, but fluctuations caused by 
a finite velocity are expected to be irrelevant, provided 
there is a surface ordering in equilibrium. The surface 
critical behavior at equilibrium is a complicated phe- 
nomenon [2514271 ] in which one should take into account 
the effect of enhanced or reduced surface couplings and 
study the critical behavior at the different fixed points 
(ordinary, extraordinary, special, surface, etc.). Also dif- 
ferent considerations should be made for systems having 
an order parameter with continuous symmetry, such as 
the XY- or the Heisenberg model [28]. Finally, one could 
also think about using different geometries, such as edges, 
wedges 0, |3(| , or parabolic shapes [3(3, El to model a 
tip sliding on a flat surface. 
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Appendix A: Self-consistency equation in ID 



Thus, we have 



Here we consider the ID Potts model in Eq. (|2.2j) in 
the presence of a field 



N 



N 



pn [s) = -kJ2s( Si - s i+ i)- h f J2^) (Al) 



and use the periodic boundary condition, s^v+i = s i- 
The transfer matrix of the problem is a q x q symmetric 
matrix (q > 2) 



(A2) 



in terms of which the partition function reads Z eq = 
TrjT^j. In the large- N limit, Z cq = A^, where A m 
is the leading eigenvalue of the transfer matrix given by: 
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(A4) 

and the magnetization follows from Eq. (I4.14[) . 

Now taking the value of the static effective field in 
Eq. (|4.12l) and using (fMJ) in Eq. (|4.14[) we obtain 



q{e K fl 



2) 



2(q-l)m f . cq -q+2 = — -= 

^ {e K n-q+2Y + A{q-l)(n + l) 

(A5) 

where f2 = S!(m/, Kf) is defined in Eq. (|4.12[) . To obtain 
a self-consistent solution, m./, C q = m fi we first take the 
square of both sides of Eq. (| A5|) leading to 



[(e K n 



-4(g-l)(fi + l)] 



x [2(q-l)m f -q+2] =q 2 (e K n-q + 2) 2 , (A6) 

which has always the trivial solution m/ = 0. Then we 
multiply Eq. (IA6|) by (1 — m/r) 2 /qrrif and obtain the 
cubic polynomial which is given in Eq. (|4. 18[) . Note that 
this polynomial has an extra root which is a solution of 
the squared equation, but not of the original one. 



Appendix B: Surface magnetization of the quantum 
Potts model for large-g values 

In the following we consider the ground state of the 
quantum Potts model defined in Eq. (I4.30[) , extended by 
the surface field term in Eq. (I4.31[) . We use fixed-spin 
boundary conditions at j = L and the surface at j = 1 is 
free. In the large-g limit, at the critical point J c = 1, the 
ground state of the system is i-fold degenerate, having 
an energy En = —J(L — 1). In a first-order perturbative 
treatment [l2T[l4| . corrections of the order of 1/y/q are 
obtained through the solution of the following secular 
eigenvalue problem, hv Q = e Q v Q . Here h is a symmetric 
L x L matrix 



(h 1 
1 8 



h = 



1 



\ 



9 1 
1 26» 



(Bl) 



1 (L-l)9j 



where plays the role of a reduced temperature, since 
the coupling is parametrized as J = 1 — 9/y/q. Using the 
components of the ground-state eigenvector vn(k) with 
k = 1, . . . , L one can express the magnetization profile 
as 



= EMfc)] 2 , 



(B2) 



fc=i 



and the surface magnetization toi = to/. 
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The surface magnetization as a function of the reduced 
temperature 9 is shown in Fig. |9] for different values of 
the surface field h. 

At the critical point, 9 — 0, the surface magnetization 
is vanishing for h < 1 and it starts linearly for small 9 
as: 



m 1 (9,h) = _ n2 +0[9 2 ae)}, h<l. (B3) 



The second derivative, d 2 mi(9, h)/d9 2 ~ diverges 
as £ ~ For ft, > 1, there is a finite surface magne- 

tization at the critical point and a small-0 behavior given 
by: 



FIG. 9: (Color online) Surface magnetization of the large-g 
state Potts model as a function of the reduced temperature 
8 for different values of the surface field h. In the disordered 
phase, 9 > 0, the surface magnetization vanishes for any finite 
value of h. For the different curves h increases from left to 
right. 



mi (9, h) = 1 - — + _ (-9) + 0(9 2 ), h>l. 

(B4) 

In the disordered phase, 9 > 0, the surface magnetization 
vanishes for any finite value of h. Consequently, there is 
a first-order surface phase transition at h = 1. 
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